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Rotation of a body, according to Einstein's theory of general rel- 
ativity, generates a "force" on other matter; in Newton's gravita- 
tional theory only the mass of a body produces a force. This phe- 
nomenon, due to currents of mass, is known as gravitomagnetism 
owing to its formal analogies with magnetism due to currents of 
electric charge. Therefore, according to general relativity, Earth's 
rotation should influence the motion of its orbiting satellites. In- 
deed, we analysed the laser ranging observations of the orbits of 
the satellites LAGEOS and LAGEOS II, using a program devel- 
oped at NASA/GSFC, and obtained the first direct measurement 
of the gravitomagnetic orbital perturbation due to the Earth's ro- 
tation, known as the Lense-Thirring effect. The accuracy of our 
measurement is about 25 %. 

In general relativity 1 ' 2 the concept of inertial frame has only a local mean- 
ing and a local inertial frame is "rotationally dragged" by mass-energy cur- 
rents, in other words moving masses influence and change the orientation of 
the axes of a local inertial frame (gyroscopes). Thus, an external current 
of mass, such as the spinning Earth, "drags" and changes the orientation of 
gyroscopes. This is the "rotational dragging of inertial frames", or "frame- 
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dragging" ("dragging", as Einstein named it). The NASA Gravity Probe-B 
experiment 3 is aimed to measure with great accuracy this phenomenon on 
the orientation of the axis of spin of a small orbiting gyroscope. However, 
the whole orbital plane of a satellite is itself a kind of enormous gyroscope 
dragged by the gravitomagnetic field (Figs. 1 and 2). Indeed, in addition to 
the rotational dragging and precession of a test gyroscope due to the angular 
momentum J of a central object, the orbit of a test particle around a central 
body with angular momentum J has a secular rate of change of the longi- 
tude of the line of the nodes (intersection between the orbital plane of the 
test particle and the equatorial plane of the central object), discovered by 
Lense-Thirring (1918): (f ense ~ imrrm a = 2J/[a 3 (1 - e 2 ) 3 / 2 ]; where a is the 
semimajor axis of the test particle, and e its orbital eccentricity. The orbit 
of the test particle also has a secular rate of change of the mean longitude of 
the orbit and of the longitude of the pericenter cD, (defining the Runge-Lenz 
vector): ~ Lense - Thm 9 = 2 J ( j _ 3 cos I l)/[a 3 (l - e 2 ) 3 / 2 ]; where 1 is the 
orbital angular momentum, unit vector, of the test particle, and / its orbital 
inclination (angle between the orbital plane and the equatorial plane of the 
central object). 

Since 1896 several experiments have been discussed and proposed to mea- 
sure the rotational dragging of inertial frames by a spinning body 1-4 . So far, 
the only indirect astrophysical evidence for the rotational dragging of inertial 
frames by a current of mass was given by the periastron precession rate of 
the binary pulsar PSR 1913+16 5 . 

Our direct measurement of the Lense-Thirring effect was obtained by 
laser ranging observations of the satellites LAGEOS and LAGEOS II 6 (Fig. 
3). The gravitomagnetic field has changed the point of closest approach to 
Earth - perigee - of the satellite LAGEOS II by about 11 meters during our 
period of observation of about 3.1 years. The semimajor axis of LAGEOS 
is a = 12, 270 km, the period P = 3.758 hr, the eccentricity e = 0.004, 
and the inclination / = 109.9°. The semimajor axis of LAGEOS II is an = 
12, 163 km, the eccentricity ejj = 0.014, and the inclination /// = 52.65°. 

Cornerstones of our analysis were: the NASA launch in 1976 of the LA- 
GEOS satellite; the development by the NASA/Goddard Space Flight Center 
of the powerful program GEODYN, and of its new version GEODYN II, for 
satellite orbit determination, geodetic parameter estimation, tracking instru- 
ment calibration, satellite orbit prediction and other applications in geodesy; 
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the determination of highly accurate Earth's gravity field solutions, includ- 
ing GEML1, GEML2, GEMT1, GEMT2, GEMT3, GEMT3S. The other new 
basic elements that made our direct measurement of the Lense-Thirring ef- 
fect possible were: (1) the launch in October 1992 by NASA and ASI of the 
laser-ranged satellite LAGEOS II; (2) the new Earth's gravity field solutions 7 
JGM-2 and JGM-3, jointly developed by NASA-Goddard and by the CSR 
(Center for Space Research) of the University of Texas at Austin; (3) the 
continuous laser ranging to the satellites LAGEOS and LAGEOS II from 
several stations around the world, the ranging data from the best stations 
have a precision of a few millimeters; and (4) the use of a new method 8 to 
measure the gravitomagnetic field. 

We analysed the orbits of the satellites LAGEOS and LAGEOS II us- 
ing existing laser ranging observations, a highly accurate modeling of their 
orbital perturbations including the gravity field solution JGM-3, and the 
1994 version of GEODYN II. All the general relativistic perturbations due 
to the masses of Earth and Sun, including the de Sitter or geodetic effect 
(today measured with accuracy of the order of 10~ 2 ), were incorporated in 
the GEODYN equations of motion and then computer-integrated; we did not 
however include in our model the orbital perturbations due to the Earth's 
angular momentum, that is the Lense-Thirring, gravitomagnetic, effect to be 
determined. In order to measure the frame-dragging effect from our residu- 
als we introduced a new parameter /i, which, by definition, is one in general 
relativity, ji GR = 1, and zero in Newtonian theory 2 . 

The residuals of the orbital elements of a satellite give a measure of any 
perturbation that is not modeled accurately enough or that is not included 
in the model. The orbital elements we analysed are: the node of LAGEOS 
I, the node of LAGEOS II, and the perigee of LAGEOS II. The nodes of 
LAGEOS and LAGEOS II are both dragged by the Earth's angular momen- 
tum; according to the Lense-Thirring formula one has: Q^ ense - Thirrm 9 ^ 

31 milliarcsec/yr and Qj T 9 = 31.5 milliarcsec/yr. The argument 

of pericenter (perigee in our analysis), u, of a test particle, that is the an- 
gle on its orbital plane measuring the departure of the pericenter from the 
equatorial plane of the central body, also has a Lense-Thirring drag; for LA- 
GEOS I one has: ^ense-T^mng ^ 32 milliarcsec/yr, and for LAGEOS 
II: ^Lense-Thimng ^ _ m nii arcsec j y T The nodal precessions of LA- 
GEOS and LAGEOS II can be determined with an accuracy of the order of 
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1 milliarcsec/yr, or less. In fact, we obtained a root mean square of the node 
residuals of about 2 milliarcsec for LAGEOS and of about 3 milliarcsec for 
LAGEOS II, over a total period of observation of about 3.1 years. Regarding 
the perigee, the observable quantity is eau, where e is the orbital eccentricity 
of the satellite. Thus, for LAGEOS the perigee precession Co is an extremely 
difficult quantity to measure; its orbital eccentricity is in fact about 4xl0~ 3 . 
The orbit of LAGEOS II is more eccentric: its orbital eccentricity is about 
0.014, and the Lense-Thirring drag of the perigee of LAGEOS II is almost 
twice as large, in magnitude, as that of LAGEOS. In fact, we obtained a 
root mean square of the residuals of the LAGEOS II perigee of about 35 
milliarcsec over about 3.1 years, whereas the total effect of frame-dragging 
on the perigee, over about 3.1 years, is = —176 milliarcsec. 

The most critical source of error in our measurement arises from uncer- 
tainties in the Earth's even zonal harmonics and in their temporal variations. 
Using only the satellites orbiting today, one cannot eliminate the unmodeled 
orbital perturbations due to all the even zonal harmonics; in particular, un- 
modeled orbital effects due to the harmonics of lower order are of a size com- 
parable to or larger than the Lense-Thirring effect. However, by analysing the 
JGM-3 solution with its uncertainties in the even zonal harmonic coefficients, 
and by calculating the secular effects of these uncertainties on the orbital ele- 
ments of LAGEOS and LAGEOS II, we found that the main sources of error 
in the determination the frame-dragging effect are concentrated in the first 2 
or 3 even zonal harmonics, that is J 2 , J4 and J 6 . To further test 10 the order 
of magnitude of the real errors in the estimated value of the J 2n coefficients, 
we took the difference 8 between two different gravity field solutions: JGM-3 
and GEMT-3S 9 . We then found that by far the largest uncertainties, on the 
nodes of LAGEOS and LAGEOS II and on the perigee of LAGEOS II, arise 
from AJ 2 and AJ 4 , a smaller error is due to AJ 6 , and much smaller errors 
arise from the differences in the other J 2n coefficients. However, we have the 
three observable quantities: the node of LAGEOS, the node of LAGEOS II, 
and the perigee of LAGEOS II, and we want to determine the parameter /i, 
measuring the frame-dragging effect. Then, we can use these three observable 
quantities Qi, flu and un to determine /i thereby eliminating the two largest 
sources of error arising from the uncertainties in J 2 and J 4 . This new method 
leads to a value of fi unaffected by the errors due to <5J 2 and <5J 4 , by far the 
largest, but sensitive only to the smaller errors due to <5J 2n with 2n > 6. As 
regards tidal, secular and seasonal changes in the geopotential coefficients, 
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we stress that the main effects on the nodes and perigee of LAGEOS and 
LAGEOS II due to tidal and other temporal variations in the Earth's gravity 
field are due to changes in the first two even zonal harmonic coefficients, J 2 
and J 4 . Any tidal error in J 2 and J 4 , and any error due to other unmodeled 
temporal variations in J 2 and J 4 , including their secular and seasonal varia- 
tions, is eliminated using our combination of residuals of nodes and perigee. 
In particular, most of the errors due to the 18.6 year and 9.3 year tides, asso- 
ciated with the Moon node, are eliminated in our measurement. Thus, using 
three observable quantities, the two nodes and the perigee, one can solve 
for fi and eliminate 5J 2 and 5J 4 : 5 tl E L ^ eosl + h5 (lT a p geosII + k 2 S ^Lageosii = 
fx ( 31 + 31.5 k\ — 57 k 2 ) milliarcsec/yr + [contributions from 5Jq, 5J%, ... ], 
where k\ = 0.295 and k 2 = —0.35 are obtained (in order to eliminate the SJ 2 
and SJ4 errors) from the system of the three equations for the nodal rates 
of LAGEOS and LAGEOS II and for the perigee rate of LAGEOS II. The 
best fit lines of the residuals of the nodes of LAGEOS and LAGEOS II had 
a slope of respectively = — 11 milliarcsec and = 40 milliarcsec, and the best 
fit line of the residuals of the perigee of LAGEOS II had a slope of = — 
188 milliarcsec. In Fig. 4 we plotted the sum of the residuals of the nodes 
of LAGEOS and LAGEOS II and perigee of LAGEOS II according to our 
formula to eliminate the 5J 2 and <5J 4 errors and after having removed 10 
small periodical residual signals (corresponding to 9 main tidal effects and 
to the largest solar radiation pressure perturbation) and the small observed 
inclination residuals. In other words each point of Fig. 4 was obtained by 
one residual of the node of LAGEOS, plus the corresponding residual of the 
node of LAGEOS II times the factor 0.295, plus the corresponding residual 
of the perigee of LAGEOS II times the factor —0.35. By fitting a straight 
line through these combined residuals of nodes and perigee (obtained using 
the JGM-3 gravity field model) we finally found: 

^ = 1.1, (1) 

This combined, measured, gravitomagnetic perturbation of the satellites' 
orbits corresponds to about 12 meters at the LAGEOS altitude, that is about 
205 milliarcsec. The root mean square of the post-fit combined residuals is 
about 13 milliarcsec. The main error sources affecting the nodes of LAGEOS 
and LAGEOS II and the perigee of LAGEOS II are: errors due to uncertain- 
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ties in the even zonal harmonics, J2„ (with In > 6 in our measurement); 
errors due to unmodeled tidal perturbations and other temporal variations 
in the Earth's gravity field (due to J 2n with 2n > 6 in our measurement), 
random and stochastic observational errors; errors due to uncertainties in 
the orbital inclinations (though we corrected nodes and perigee with the 
residuals of the orbital inclinations); errors due to nongravitational pertur- 
bations, including direct solar radiation pressure, Earth's albedo, Yarkovsky 
anisotropic thermal radiation, Rubincam effect (anisotropic re-radiation of 
Earth infrared radiation absorbed by the LAGEOS retro- reflectors) , particle 
drag, and errors due to the estimated values of the satellite's reflectivities 
and estimated 15-day along track accelerations. 

By calculating the effects of all these systematic and random error sources 
(paper in preparation), we found: 

5^ 25% a (2) 

In conclusion, we obtained the result: 

u = 1.1 ±0.25 (3) 
(whereas a — 1 in general relativity). 
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FIG. 1 The gravitomagnetic field, H, generated by the angular momentum, 
J, of a central rotating body. In general relativity, for a localized, stationary, 
mass-energy distribution, in the weak- field and slow- motion limit, the three 
"vector" components of the metric tensor are given by: h = — 2(3 x x)/r 3 , 
where J is the angular momentum of the central body and h is known 
as the gravitomagnetic potential. The gravitomagnetic field H is given by 
H = V x h. To characterize the gravitomagnetic field generated by the an- 
gular momentum of a body and the Lense-Thirring effect, and distinguish 
it from other relativistic phenomena such as the de Sitter effect - due to 
the motion of a gyroscope in a static gravitational field - one may give a 
description of the gravitomagnetic field in terms of spacetime-curvature in- 
variants. The pseudoinvariant *R • R, built from the Riemann tensor R and 
its dual *R, gives an invariant characterization of gravitomagnetism since it 
is nonzero in the field of a central body if and only if the body is rotating. 
Indeed the pseudoinvariant *R- R is proportional to the angular momentum 
of the central body. Thus, one may describe gravitomagnetism as that phe- 
nomenon of nature such that spacetime curvature is generated by the spin of 
a body 2 . 

FIG. 2 A twisted jet from the nucleus of the galaxy 3C 66B. This ultraviolet 
picture has been taken by the ESA Faint Object Camera of the NASA Hub- 
ble Space Telescope. The twisted jet of plasma extends 10,000 light-years 
from the nucleus of the galaxy 3C 66B located at about 270 million light- 
years from Earth. The ultraviolet radiation is emitted by electrons in the 
jet spiraling through magnetic fields. Long radio jets from quasars and ac- 
tive galactic nuclei are observed to have constant directions in space, which 
correspond to emission time scales that may reach millions of years. The 
constant direction of the jets suggests the existence of a central astrophysi- 
cal gyroscope; this engine and gyroscope might be a super-massive spinning 
black hole with its gravitomagnetic field. The constant orientation of the 
emitted jets may then be explained 11 using the gravitomagnetic field of the 
central spinning body. 

FIG. 3 The laser-ranged satellite LAGEOS II. Laser ranging to the Moon and 
to artificial satellites is an impressive technique to measure distances from 
a laser-tracking station on Earth to retro-reflectors placed on the Moon, or 
on satellites orbiting Earth. By the use of short laser pulses ranges can be 
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measured with accuracies of less than 1 cm from emitting lasers on Earth 
to retro-reflectors on a satellite, and with accuracies of less than 10 cm to 
retro-reflectors on the Moon. The NASA-ASI (Italian Space Agency) satel- 
lite LAGEOS II is a high- altitude, small cross-sectional area-to-mass ratio, 
spherical, laser-ranged satellite. It is made of heavy brass and aluminum, is 
completely passive and covered with laser retro-reflectors. It acts as a ref- 
erence target for ground-based laser-tracking systems to measure - via laser 
ranging - crustal movements, plate motion, polar motion and Earth rotation. 
LAGEOS II is essentially identical to the NASA satellite LAGEOS (LAser 
GEOdynamics Satellite) but they have different orbital parameters. 

FIG. 4 Sum of the residuals of the nodes of LAGEOS and LAGEOS II 
and perigee of LAGEOS II from November 1992 to December 1995, us- 
ing the method described in the text. On the vertical axis we plotted 
[node residuals of LAGEOS) + 0.295 (node residuals of LAGEOSII) — 
— 0.35 (perigee residuals of LAGEOSII). In our analysis we included polar 
motion from VLBI (IERS), Earth's solid and ocean tides and Earth's gravity 
field, GM and spherical harmonics up to order 50, from the JGM-3 gravity 
field model, solar, lunar and planetary perturbations and nongravitational 
perturbations including solar radiation pressure, Earth's radiation pressure, 
anisotropic thermal radiation effects, and atmospheric drag. For each 15-day 
arc we estimated all station coordinates except the latitude of Goddard Space 
Flight Center and the latitude /longitude of Hawaii (maintained fixed), the 
spacecrafts' initial conditions (initial positions and velocities), the satellites' 
reflectivities and 15-day along-track accelerations. The best fit line shown 
through these combined residuals has a slope of about 66 milliarcsec/yr (the 
total integrated effect corresponds to about 12 meters at the LAGEOS alti- 
tude), that is ji = 1.1 (whereas /i = 1 in general relativity), and the corre- 
sponding root mean square of the residuals is about 13 milliarcsec. Due to 
systematic (secular and periodical) errors and random errors, we estimated 
the total error in our measurement of ji to be less than 25% of \x. 
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Figure 4 
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